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^q", Abstract 

' Wc present a nonparametric method to determine the sign of w + 1 in the equation of state of 

{SJ . dark energy. It is based on geometrical behaviour and is more tolerant to uncertainties of other 

' cosmological parameters than fitting methods. It permits to distinguish between different classes 

. of dark energy models even with relatively low precision data. We apply this method to SNLS 

Q ' supernovae and to gold sample of re-analyzed supernovae data from Riess et a/.2004 [Tj. Both data 

sets show strong indication of w < — 1. If this result is confirmed by more extended and precise 
data available in near future, many of dark energy models, including simple cosmological constant, 
, standard quintessence models without interaction between quintessence scalar field(s) and matter, 

J> ' and scaling models are ruled out. 

o : 

in . 

, Recent observations of Super Novae (SN), Cosmic Microwave Background (CMB), and Large Scale 
Structures (LSS) indicate that the dominant content of the Universe is a mysterious energy with an 
equation of state very close to Einstein Cosmological Constant. The equation of state is defined by w, 
the ratio of pressure p to density p, w = Pj p. For a cosmological constant w = —1. The observed mean 

\ value of w for dark energy is very close to —1. Some of the most recent estimations are the followings: 
Qh 1 From combination of 3-year WMAP and SuperNova Legacy Survey (SNLS), w = — 0.97(]~q 9 07 [2]; from 
q . combination of 3-year WMAP, large scale structure and supernova data, w = -1-061^09 [2]; from 

| combination of CMAGIC supernovae analysis and baryon acoustic peak in SDSS galaxy clustering 
statistics at z = 0.35, w = — 1.21 jl ; 12 [3]; and finally from baryon acoustic peak alone w = — 0.8±0.18. 
It is evident that with inclusion of one or two sigma uncertainty to measured mean values, the range 
of possible values for w runs across the critical value of —1. Moreover, in all these measurements 
the value of w depends on other cosmological parameters and their uncertainties in a complex way. 
Reconstruction methods for determining cosmological parameters from observations (see [4] and 
references therein for a review of methods) usually use fitting of continuous parameters on the data 
and determine a range for w. 

On the other hand, the sign of 7 = w + 1 is more crucial for distinguishing between various dark 
energy models than its exact value. For instance, if 7 < 0, scalar field (quintessence) models with 
conventional kinetic energy and potential are ruled out because in these models 7 is always positive. 
Decay of dark matter to dark energy [7] [8], or in general an interaction between these components 
can lead to an effective 7 < without violating null energy condition [9] [lOj . 

Here we propose a nonparametric method specially suitable for estimating the sign of 7. When 
the quality of data is adequate, the quantity A(z) defined in ([2]) can also be used to fit the data 
and to measure the value of 7. The expression nonparametric here is borrowed from signal processing 
literature where it means testing a null hypothesis against an alternative hypothesis by using a discrete 
condition such as jump, sign changing, etc., in contrast to constraining a continuous parameter (see 
e.g. [II]). We show that geometrical properties of A(z) are related to the sign of 7 and we can detect 
it without fitting a continuous parameter. 



The density of the Universe at redshift z is: 

n m (i + z f + n h (i + zf + n de (i + z 



where p(z) and pq are total density at redshift z and in local Universe, respectively; Q m , Q^, and f^e are 
respectively cold and hot matter, and dark energy fraction in the total density at z = 0. We consider a 
flat universe in accordance with recent observations [2]. At low redshifts, the contribution of CMB to 
the total mass of the Universe is negligible. The contribution of neutrinos is £l v h 2 = eV, 
h = H^/lOOkm Mpc~ 1 sec~ 1 , where Hq is present Hubble constant. The upper limit on the sum of 
masses of neutrinos from 3-year WMAP is ^2m u < 0.62 (95% confidence level) [13]. Therefore, for 
z < 1 their contribution to the total mass of the Universe is < 4% even if one of the neutrinos has 
very small mass and behaves as a warm dark matter. This is less than the uncertainty on the fraction 
of dark matter, and thus the approximation Q m + f^ e ~ 1 is justified. It can be easily shown that in 
this case: 

3(1 + z) z po dz 

Similar expressions can be obtained for non-standard cosmologies such as DGP |14j model and other 
string/brane inspired cosmologies [15]. It is also possible to find an expression similar to ([2]) for non- 
flat FLRW models and without neglecting hot matter. The left hand side would however depend on 
fife, and £lh and would be more complex. Nonetheless, when the contribution of these components at 
low redshifts are much smaller than cold matter and dark energy, the general behaviour of A{z) will 
be the same as approximate case studied here. 

The right hand side of ([2]) has the same sign as 7. Moreover, the sign of its derivative is opposite to 
the sign of 7 because due to the smallness of observed 7, the term 7 — 1 is negative. This means that 
A(z) is a concave or convex function of redshift, respectively for positive or negative 7, see figure [TJ 
In the case of a cosmological constant A(z) = for all redshifts. This second feature of expression 
([2]) is interesting because if fi TO is not correctly estimated, A(z) will be shifted by a constant, but this 
will not modify geometrical properties of A{z). 

The left side of expression ([2]) can be directly estimated from observations. More specifically, f2 m is 
determined from conjunction of CMB, LSS, and supernova type la data, and at present it is believed 
to be known with a precision of ~ 5%. At low redshifts, the derivative of the density is best estimated 
from SN type la observations. In the case of FLRW cosmologies, the density and its derivative can be 
related to luminosity distance Di, and its first and second derivatives: 

1 A 2 (dDi _ _Dj^\ _ d 2 Di 

«/ A _ 1 a P _ 1 + Z\ dz 1 + Z> ~dW /o\ 

"W- 3(1 + 1(^-^) 3 [ ) 

Di = (1 + z)H [ , H\z) = ^p(z) (4) 

It is remarkable that the right hand side of ([3]) depends only on one cosmological parameter, Hq. 
Nonetheless, similar to an uncertainty on Q m , Hq scales B{z) similarly at all redshifts, and therefore 
it does not change the overall geometrical properties B(z) and A(z). D\ can be directly obtained from 
observed luminosity of standard candles such as supernovae type la. In the case of LSS observations 
where the measured quantity is the evolution of density p with redshift, B(z) is measured directly up 
to an overall scaling by po. This does not change the geometrical properties of B(z) and A(z), i.e. the 
scaling by a positive constant does not flip a convex curve to concave or vis-versa. 

In summary, uncertainties of Q m and Hq do not affect the detection of the sign of 7 through geometrical 
properties of A(z). This is quite different from fitting methods. They are sensitive to all numerical 
parameters Hq, O m , Orf e , and id in a complex way, usually through a non- linear equation such as 
chi-square or likelihood equation, and it is very difficult to assess the effect of uncertainty of one 
parameter on the estimation of others, and more specifically on the determination of the sign of 7. 

When this sign detection method is applied to standard candle data such as supernovae type la where 
the measured values are related to Di , according to ([3]) , one has to calculate first and second derivative 
of Di. Numerical calculation of derivatives is not trivial. To have a stable and enough precise result, 
not only the data must have high resolution and low scatter, but also it is necessary to smooth them. 



To test the stability of numerical calculation and the method in general, we also apply it to simulated 
data. The details of numerical methods is discussed in the Appendix. 
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Figure 1: A{z) in Eq.(|2|) as a function of redshift; a: 7 = ±0.06 (corresponding to best fit of 3-year 
WMAP, LSS, and SN data [2]) and b: 7 = ±0.2 (corresponding to best fit of 3-year WMAP and 
CMAGIC SN data [3]). For both plots H = 73 km Mpc' 1 sec' 1 and Q de = 0.77. We consider 5% of 
uncertainty for both quantities. Blue curves present right hand side of ([2]), including the uncertainty on 
the value of tide- Curves decreasing or increasing with increasing redshift present respectively positive 
and negative 7. Magenta and purple curves present numerical determination of A(z) respectively for 
positive and negative 7 from about 150 simulated supernovae distributed with uniform probability in 
log(z). These plots show that numerical errors dominate near two boundaries of the redshift range. 
It is also evident that with this method one can distinguish between positive and negative 7 for 



Two examples of reconstruction of simulated data are shown in FigHJ It is evident that despite 
deformation of the reconstructed curves due to numerical errors and uncertainties, the difference 
between convexity of curves for positive and negative 7 is mostly preserved and can be used visually 
or by using a slope detection algorithm to find the sign of 7. The simulated data is however much more 
uniform than presently available data, see Figj2j Although for mid-range redshifts the data follow a 
curve similar to what is expected from FLRW cosmologies, artifacts appear close to boundaries and 
at high redshifts where the quality of data is worse. Moreover, visual inspections or slope detection 
lack a quantitative estimation of uncertainty of measured sign for 7. Another complexity of this 
cosmological sign detection problem is that not only the observable A{z) is noisy, it also varies. In 
signal processing, in most practically interesting cases the signal is constant but noisy. Therefore, 
usual binomial estimation of the probability or optimization of detection [12] are not applicable. 

Here we take another strategy, specially suitable for this cosmological sign detection task. The null 
hypothesis for dark energy is 7 = 0. Assuming a Gaussian distribution for uncertainty of reconstructed 
A(z) from data and from simulated data for 7 = model, for each data-point we calculate the 
probability that the data-point belongs to the null hypothesis. To include the uncertainty of data, we 
integrate the uncertainty distribution from —a to +a around the mean value: 



where A% and Oi are from i data-point, and A and uqi from simulated null hypothesis model at the 
same redshift. Averaging over gives P, an overall probability that the dataset corresponds to the 
null hypothesis. As 7 = is the limit case for 7 > 0, P is also the maximum probability of 7 > 0. 

We have applied this sign detection algorithm to two supernova datasets: published data from Super- 
nova Legacy Survey (SNLS) [5] and low redshift supernovae (z < 0.45) of gold sample of re-analyzed 



M > 0-06. 





supernovae data by Riess et a/.2004 pQ. The reason for using only the low redshift subset of the latter 
compilation is that the scatter and uncertainty of the peak magnitude at higher redshifts is too large, 
and with numerical methods used in this work, it is not possible to recover a reasonably stable and 
smooth distribution for A(z). 

Figj2] shows A{z) obtained from these data. To estimate the effect of the reconstruction, we compare 
A(z) from data with simulated data as described in the caption of Fig[2j Simulated standard sources 
are at the same redshifts as in the datasets to make simulated samples as similar to data as possible. 
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Figure 2: A(z) from: 117 supernovae of the SNLS data (a-purple curve), and 88 supernovae with 
z < 0.45 recompiled and re-analyzed by Riess et aZ.2004 pQ (b-purple curve). Error bars present 
1-sigma uncertainty. In (a), green, orange, yellow, and light green curves present the reconstruction 
of A(z) from simulations for 7 = —0.2, —0.06, 0.6, 0.2, respectively. The overall probability of null 
hypothesis (7 = 0) is P = 0.27, therefore the probability of 7 < 0, 1 — P = 0.73. Light grey and cyan 
curves are theoretical calculation including the uncertainty on f^ e , respectively for 7 = ±0.06, ±0.2. 
In (b), the pink curve presents simulated distribution for 7 = - a cosmological constant. For this 
dataset 1 — P = 0.75. The dark grey straight line is the theoretical expectation for A(z) when dark 
energy is a cosmological constant. Definition of other curves are the same as (a). 

In both datasets the probability of 7 < or equivalently w < — 1 is larger than 70%. The SNLS data 
is consistent with a 7 as small as ~ -0.50. There are however significant deviations from a smooth 
distribution for z < 0.1 and z > 0.5. We attribute them to relatively large scatter of the data at these 
redshifts that makes reconstruction instable, see Figj3}-a. 

To see whether the large negative 7 concluded from the SNLS data is due to the data scattering and /or 
reconstruction algorithm artifacts, we also apply the same formalism to a subset of these data with 
z < 0.45. The result is shown in Figj3]-b along with simulations in the same redshift range. The bump 
at very low redshift in Figl2]does not exist in this plot, and therefore we conclude that it is an artifact 
of numerical calculation. Although A(z) distribution in this data set is also convex and the probability 
of 7 < is > 90%, it does not have the same slope as any of models. More specifically, it seems that 
low and high redshift sections of the curve correspond to different values of 7. For z < 0.15, A(z) is 
close to theoretical and simulated data with 7 = —0.2. For z > 0.25, A(z) approaches the values for 
larger and even positive 7. Such behaviour does not appear in Figj2j Giving the fact that the number 
of available data points with z > 0.25 in this subset is small the most plausible explanation is simply 
numerical artifacts. Alternative explanations are evolution of 7 with redshift and an under-estimation 
of Q m , see Eq.([2]). If the latter case is ture, the value of 7 must be even smaller than —0.2. With a data 
gap in 0.15 < z < 0.25 interval and a small total number of entries in this data set- only 58 supernovae 
- it is not possible to make any definite conclusion about the behaviour of this data. We should also 
mention that nearby SNe with z < 0.25, intermediate SNe 0.25 < z < 0.4, and high redshifts ones 



*As the main purpose of this paper is determination of the sign of 7, we don't perform any fitting to obtain its values. 
Estimation of values here are from plots. 



z > 0.4 in SNLS are not completely treated in the same way [5]. It is therefore possible that some of 
the observed features explained here are purely artifacts of the analysis of raw observations. 




+ 



iiiiii i i 

i! 



inn 1 1 



0.2 0.3 
rcdahift 



Figure 3: a: /i;, distribution of SNLS supernovae, data (magenta), smoothed distribution by the 
method explained in the Appendix (blue). Although this distribution look quite smooth, even small 
sudden variations can lead to large variations in derivatives, b: A(z) for SNLS supernovae with 
z < 0.45. definition of curves is the same as in Fig|2j For this dataset 1 — P = 0.93 when cosmological 
parameters are the same as Figj2j If = 0.73 is used, 1 — P = 0.96 for the same data. 

Fig|2}b shows A(z) determined from gold sample supernovae recompiled by Riess, et al. pp. This 
dataset is consistent with 7 < with a probability ~ 75% at 1-sigma and ~ 66% at 2-sigma. This 
plot also shows that the value of 7 estimated from this data is ~ —0.06, larger than estimation from 
SNLS data. 

The reason for the difference between estimated values for 7 from SNLS and Riess, et decompilation 
is not clear because the value of w reported in Ref. [1] and Ref [5] are consistent. Nonetheless, the 
estimated values here are in the range reported by other works [2] [3]. The difference between results 
of two datasets is most probably related to their different scatter and uncertainty, and the fact that 
low and high redshift data are not treated in the same way. We should also mention that recent claims 
about contamination of supernova type Ic and the effect of asymmetric explosion in the lightcurve 
of supernova type la [16], and possible differences between low and high redshift supernovae [17] can 
not explain the difference between the results of these datasets. These phenomena affect both surveys 
in the same way. Despite these discrepancies and uncertainties, both datasets are in good agreement 
about negative sign of 7. 

The results of present study show various shortcomings in both datasets used here. Our first remark 
is the large gap in redshift distribution of observed supernovae in redshift range 0.1 < z < 0.3. Both 
datasets have less than 6 supernovae in this range and nothing in 0.15 < z < 0.25. This is not an 
important issue for finding redshift distribution of Di, but redshift gap becomes very important when 
derivatives of Di are calculated. The lack or rareness of supernovae data in this redshift range is related 
partly to the history of star formation in galaxies [18] , and partly to optimization of surveys |19j for 
detection of very low or high redshift supernovae which decreases the probability of detection of mid- 
range SNe. A systematic survey of galaxies in this range is therefore necessary to fill the present gap. 
Sloan Supernova Survey [20] is optimized to detect SNe in 0.1 < z < 0.35 and should provide the 
missing data in near future. 

Our second observation is a large scatter in both datasets around redshift ~ 0.55, see Fig|3]-a. This 
leads to a large scatter in the numerical determination of derivatives in ([3]) and makes the results 
unusable, see Fig|2}-a. In future observations the reason of this large scatter should be understood, 
and if possible reduced. From theoretical calculation and simulations shown in Figj2] one can also 
conclude that with present uncertainties of cosmological parameters, the most important redshift range 
for determining the equation of state of dark energy is z < 0.8. Although the technical challenge, 



understanding of physics of supernovae [T7] and their evolution, and applications for other astronomical 
ends make the search of supernovae at larger distances interesting, they will not be in much use for 
determining the equation of state of dark energy, at least not at lowest level which is the determination 
of redshift independent component of w. 

On the other hand, improvements in numerical techniques and algorithms would lead to better mea- 
surements. One of the possibilities in this direction is the application of an adaptive smoothing 
algorithm with variable degrees of smoothing depending on the amount of scattering in the data. 
More sophisticated smoothing algorithm also has been proposed [21| . We postpone the application 
these more advanced methods to future, when a larger data set becomes available. 

In summary, we have proposed a nonparametric formalism to investigate the sign of 7 in the evolution 
equation of dark energy. When better quality data become available fitting can be added to to this 
method to find the value of 7 and not just its sign. The advantage of this method with respect to multi- 
parameter fitting is that the dependence on cosmological parameters is explicit, and therefore it is 
easier to assess the effect of their uncertainties on the measurements. This method is specially suitable 
for applying to supernova data where the standard observable - the peak magnitude - can be directly 
related to cosmological distance, and thereby to cosmological variation of total density. It can also 
be applied to data from galaxy clustering surveys which permit to determine the variation of average 
density of the Universe with redshift, but not to integrated observables such as CMB anisotropy. By 
applying this method to two of largest publicly available supernovae data sets we showed that they 
are consistent with a w < — 1. Present data is not however enough precise to permit the estimation 
of \w\ with good certainty. With on-going projects such as SNLS, Supernova Cosmology Project |22j, 
and SDSS SNe survey, and future projects such as SNAP and DUNE, enough precise datasets should 
be available in the near future. 

Appendix: To calculate A(z) in ([2]) for standard candles we must determine the luminosity distance. 
It is related to the magnitude of the standard candle: Di/Dq = KF''/ 5 , where Dq is the distance for 
which the common luminosity of standard sources is determined from theoretical models or observa- 
tions. The standard magnitude at a given distance depends on Hq and a correction must be added if a 
different Hq is used in the calculation of ^ [5]. For simulating the data D\ is calculated from ^ and 
the relation above is used to determine the corresponding /i& to which we add a random uncertainty 
with a standard deviation of 3%. From this point the same procedure is applied to both simulated 
and real data to determine A(z). 

Expression ([3]) for B(z) contains first and second derivatives of Di which must be calculated numeri- 
cally from data. It is however well known that direct determination of derivatives leads to large and 
unacceptable deviation from exact values. One of the most popular alternative methods to the direct 
calculation is fitting of a polynomial around each data point and then calculating an analytical deriva- 
tive using the polynomial approximation in place of the data. We use this approach to determine 
derivatives of fib and D[. In addition, before applying this approach, we smooth the distribution of 
magnitudes using again the same polynomial fitting algorithm. 

In FLRW cosmologies the redshift evolution of the luminosity distance is very smooth. Therefore, 
a second order polynomial for smoothing is adequate. Fitting is based on a right-left symmetric, 
least x 2 algorithm, and we have implemented Ifit function of Numerical Recipes [23] for this purpose. 
By trial and error we find that 19 data-point fitting gives the most optimal results regarding the 
number of available data points and their scatter. Close to boundaries however less data point for 
fitting is available in one side of each point, and therefore the fitting is less precise. The artifacts 
discussed above are mostly related to this imprecision of numerical calculation. In the present work 
no adaptive smoothing is applied. In addition to smoothed data and their derivatives, the function Ifit 
calculates a covariant matrix for uncertainty of parameters (derivatives). We use diagonal elements as 
1-sigma uncertainty of the smoothed data and its derivatives. The uncertainty of A(z) is calculated 
from the uncertainty of terms in ([2]) and ([3]) using error propagation relation: For f(xi, X2, ■ ■ ■), 
= Y^i a xi(9f/dxi) 2 . Smoothed terms and parameters in A(z) are considered as independent 



variables with their own uncertainty. 
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